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ABSTRACT 

We compute black hole masses and bolometric luminosities for 57 active galactic nuclei (AGN) in 
the redshift range 1.25 < z < 2.67, selected from the GOODS-South deep multi-wavelength survey 
field via their X-ray emission. We determine host galaxy morphological parameters by separating 
the galaxies from their central point sources in deep HST images, and host stellar masses and colors 
by multi-wavelength SED fitting. 90% of GOODS AGN at these redshifts have detected rest-frame 
optical nuclear point sources; bolometric luminosities range from 2 x 10 43 " 46 erg s" 1 . The black holes 
are growing at a range of accretion rates, with > 50% of the sample having L/LEdd < 0.1. 70% of host 
galaxies have stellar masses M* > 10 10 M Q , with a range of colors suggesting a complex star formation 
history. We find no evolution of AGN bolometric luminosity within the sample, and no correlation 
between AGN bolometric luminosity and host stellar mass, color or morphology. Fully half the sample 
of host galaxies is disk-dominated, with another 25% having strong disk components. Fewer than 15% 
of the systems appear to be at some stage of a major merger. These moderate-luminosity AGN hosts 
are therefore inconsistent with a dynamical history dominated by mergers strong enough to destroy 
disks, indicating minor mergers or secular processes dominate the co-evolution of galaxies and their 
central black holes at z ~ 2. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: fundamental parameters — galaxies: 
bulges — galaxies: Seyfcrt 



1. INTRODUCTION 



The blac k hole-galaxy connection is well-established at 
local fe.g-lFerrarese fc FordHl999t iGebhardt et aT1 [2000l: 
iMarconi fc Huntl 120031: IHaring fc Rixll2004D and moder- 



Treu et all 12004 iWoo et al 



120051: |Pen7 et al.l 



120061: iBennert et al.ll2010Ll2011aD redshifts. Clues to un- 
derstanding the co-evolution of galaxies and their cen- 
tral supermassive black holes come from detailed study 
of both moderate-luminosity growing black holes and the 
galaxies that host them. It is these moderate-luminosity 
AGN and host galaxies that cumulatively dominate the 
observed X-ray background and thus constitute a large 
fracti on of the underlying black hole growth in th e uni- 
verse (|Hasinger et al.ll2005l; ITreister fc Urry||201lD . 

Examination of both AGN and host galaxies re- 
quires separation of their blended signal, preferentially 
with a combination of space- and ground-based data at 
multiple wavelengths. Using deep Hubble Space Tele- 
scope (HST) data to extract nuclear point sources from 
AGN+host images requires care but provides the lever- 
age needed to determine AGN bolometric (total) lumi- 
nosities and host galaxy stellar masses and mo rphologies 
(jSimmons fc Urrvl[2008l : iSimmons et al.ll2011h . 

Until recently, such analyses for large, uniform 
AGN+host samples in the rest-frame optical were only 
possible to z < 1.3 using deep survey data taken with 
the Advanced Camera for Surveys (ACS) on HST. Re- 
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suits of multiple studies of host galaxy morphologies 
and black hole growth rates at those redshifts indicate 
a complex picture, with moderate-luminosity AGN pre- 
domi nantly powered by hi gh-mass, slow-growing black 
holes (jSimmons et al.ll201lT) in host galaxies with a range 
of morphologies jSanchez et all[200l iBundv et all feOOS: 
IGabor et al.ll2009j ). 

Mergers are often invoked to tr i gger black hole ac- 
cretion (e.g.. [Sanders et al.l [l988t ICroton et all |2006; 
iHopkins et all l2006al) . but mergers are equally com- 
mon in gala xies hosting active ly growing and inactive 
black holes (jGrogin et al.l 120051 ) . Moreover, moderate- 
luminosity AGN at z ~ 1, which are powered by high- 
mass black holes that have already undergone signifi- 
cant growth, common ly found in disk-dominated systems 
(jSimmons et al.l [20111 ). This raises the issue of whether 
black hole growth is s i milar at z ~ 2, wh ere star forma- 
tion dLillv et al.lll996j .iMadau et al.lll998ft an d AGN ac- 
tivitv (e.g- lFan et al.ll2001t iWall et al.ll2005j ) both peak, 
or whether mergers are relatively more important. 

The recently-installed Wide Field Camera-3 (WFC3) 
instrument on HST now provides high-resolution, rest- 
frame optical images of AGN and host galaxies to z ~ 2. 
With the Cosmic Assembly Near-in frared Deep Extra - 
galactic Legacy Survey (CANDELS; IGrogin et al.ll20ll 
imaging in the deep Chandra legacy fields, we can sepa- 
rate AGN and hosts at rest-frame optical and UV wave- 
lengths, enabling detailed investigations of the proper- 
ties of both moderate-luminosity growing black holes and 
their host galaxies. 



2 



Simmons et al. 



In Scction[2]wc describe our sample selection and anal- 
ysis methods. Section [3] presents the AGN host galaxy 
morphologies, stellar masses and colors, as well as es- 
timates of black hole masses and growth rates for the 
AGN in the sample. Section @] discusses these results 
in the overall context of black hole-galaxy co-evolution, 
specifically merger-driven versus secular growth scenar- 
ios. 

Throughout this paper, we adopt a ACDM cosmol- 
ogy with H a = 71 k m s" 1 Mpc _1 ,ft M = 0.27, A = 0.73 
(|Spergel et alJl200l . 

2. DATA 
2.1. Sample Selection 

The Chandra Deep Field - South is one of the best- 
observed parts of the sky, with deep multi-wavelength 
coverage from both ground- and space-based observa- 
tories. In particular, the 4 Megas econds of combi ned 
X-ray observations within the field (|Xue et al.ll2011l) al- 
low for selection of AGN at z ~ 2 that is relatively 
unbiased with respect to obscuration. Examination of 
those AGN and their hosts at rest-frame optical wave- 
lengths is possible with the multi-epoch deep observa- 
tions of the GO ODS-South field (chose n to coincide with 
the CDF-S; se e iGiavalisco et al.ll2004D within the CAN- 
DELS survey (jGrogin et al.ll201lD . We match the X-ray 
sources to WFC3 sources using a maximum-likelihood 
method (jCappelluti et al.ll2009t ICardamone et al.ll2010D 
that considers both positional overlap and uncertainties 
as well as multi-wavelength source fluxes to determine 
the most probable optical match to each X-ray source. 

We select AGN within the CTiandra-GOODS-S- 
CANDELS field according to the following criteria: 

1. Redshifts (compiled bv lXue e t al. 201l]) within the 
range 1.25 < z < 2.67, so that the observed wave- 
length ranges spanned by the F125W and F160W 
WFC3 HST filters are in the rest-frame B band; 

2. Source ty p e ide ntification as AGN within the 
IXue et all (|201lD 4 Ms X-ray catalog. Briefly, 
AGN are identified via X-ray luminosity, photon in- 
dex, X-ray-to-optical flux ratio, X-ray excess over 
measured star formation, and/or optical spectro- 
sc opic featur es. The criteria are described further 
in lXue et ahl ; 

3. Extended source detection within the WFC3 im- 
age with a signal-to-noise of at least 5 per pixel to 
ensure robustness of morphological fitting; and 

4. Extended source detections in at least five opti- 
cal and infrared bands (including WFC3 JH from 
CANDELS as well as ACS BVIz and ground-based 
Ks from GOODS; details of GOO DS observations 
are g iven in Tables 1 and 2 in IGiavalisco et al.l 
2004), to facilitate accurate determination of stel- 
lar masses and rest-frame colors. 

Of the 57 X-ray-selected AGN in the sample us- 
ing the abov e criteria, 26 of the redshift s are 
spectroscopic (iSzokolv et all [20041: iZheng et al.l 120041: 
Mignoli et all 120051: iVanzella et al.l l2008t iPopesso et al l 
2009t ISilverman et al.l I2OIOD and 31 a r e photometric 
(|Luo 1 et all 120101 : ICardamone et al.l 120101: iRaffertv et all 



l2011f h 41 are considered luminous AGN by IXue et al.l 
owing to their having absorption-corrected, rest-frame 
0.5 — 8 keV luminosities Lx > 3 x 10 42 erg s -1 . An addi- 
tional 14 sources have 10 42 < Lx < 3 x 10 42 erg s _1 , 
and the remaining two have luminosities of 9.1 and 
9.5 x 10 41 erg s -1 . The median and mean X-ray lu- 
minosities of the sample are 4.4 and 7.3 x 10 42 erg s _1 , 
respectively, and the most luminous X-ray AGN in the 
sample has Lx = 2.6 x 10 44 erg s _1 , consistent with ex- 
pected detections within the volu me of GOODS-S ba sed 
on the AGN luminosity fu nction (iCroom et alJl20"o"4T ). 

The CANDELS survey (jGrogin et al.l 1201 ID includes 
several deep and wide fields on the sky with both the 
WFC3 and ACS cameras on HST. We perform mor- 
phological fits on the publicly available co-added six- 
epoch images in both the F125W and F160W bands, 
which were drizzled to a resolution o f 0" .06 per pixel. 
For de tails of the image processing, see lKoekemoer et al.l 
(|20T1 . 

Morphological decomposition of two-dimensional 
AGN+host galaxy images requires accurate determina- 
tion of the point-spread function (PSF). We determine 
the PSF of the F125W and F160W CANDELS images 
using the PSF-modeling routines in the IRAF package 
daopliot. This method balances the noise- free attributes 
of a simulated PSF model with the detailed particulars 
of the actual PSF of multi-drizzled, multi-epoch data. 
We verified that the WFC3 PSF is very stable across the 
deep mosaic in the GOODS-S field, so we determine a 
single PSF for each band and use it in the morphological 
analysis of each of the sources. 

2.2. Morphological Fitting 

We perform two-dimensional parametric morphologi- 
cal fitting on the F125 W and F160W im ages for the 57 
sources using GALFIT (|Peng et al.ll200l . which is very 
reliable for recovering both normal galaxy morphologies 
and AGN- host decomposition (see lHaussler et al.l [20071 : 
ISimmons fclJrTvi l2TjQ8t iPierce et al.ll2010al for relevant 
simulations). 

Our detailed fitting procedure is similar to that per- 
formed in previ ous simulations and stud i es of AGN host 
morp hologies ([Simmons fc Urrvl 120081 : ISimmons etahl 
1201 If ) . We determine initial guesses for fit parameters 
such as total source magnitude, half-light radius, axis 
ratio, p osition angle and centro id position using SEx- 
tractor (jBertin fc Arnoutd [19911 . The value of the sky 
background varies slightly across the CANDELS mosaic, 
so we fix the SExtractor-determined value for each source 
and each iteration of the morphological fitting in order 
to avoid the confusion of sky background and extended 
source light that can ensue when the sky background is 
allowed to vary along with the other fit components. 

We fit each galaxy with two models i n each band: (1) a 
host galaxy modeled by a Sersic profile |S6rsic 1968]), and 
(2) a host galaxy with a central point source. We deter- 
mine the final fit parameters via an iterative process that 
initially fits only the central portion of the AGN+host, 
to find the best-fit centroid positions of each fit compo- 
nent. Those positions are then held fixed in successive 
iterations, which zoom outward to eventually encompass 
a large enough image region such that the extended por- 
tions of the central source are well sampled. Where 
necessary, we fit nearby, bright and/or extended com- 
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Fig. 1— HST WFC3 F125W + FIG0W images (left), 2- 
dimensional morphological fits (center), and residuals (image - fit; 
right) for four X-ray selected AGN+hosts in this sample (from top 
to bottom): a disk-dominated, intermediate, and strong-bulge host 
morphology, as well as a clumpy/asymmetric host galaxy. Object 
389 does not have a detected nuclear point source; the other three 
have point sources that are faint but detected at > 5cr. 

panions in order to simultaneously model their extended 
light distributions. Companion fitting (and masking of 
fainter, more compact companions further from the cen- 
tral source) ensures a more realistic fit to the central 
source. The first three fits to each source are performed 
with the use of a customized adaptive batch-fitting code. 
Each source is then followed up individually to ensure 
the true minimum reduced xt value within the parame- 
ter space is achieved. 

For each source and each fit (Sersic only and 
Sersic+Point Source), we perform one final iteration 
where every parameter (with the exception of the sky 
background) is allowed to vary. This final iteration en- 
sures that GALFIT has found the true best fit and also 
allows for an estimate of the computational uncertain- 
ties for each parameter, including the centroid positions 
of each component. For all sources, the differences be- 
tween the final fit parameters and those of the penul- 
timate iteration are well within the reported parameter 
uncertainties; thus, we are confident that we have found 
the true best-fit parameters for each source and fit type. 
Typically, an AGN+host galaxy in the sample requires 
6-8 total iterations for a Sersic+PSF fit, with 1-4 fit- 
ted companions. Some of the sources require additional 
iterations, usually due to the presence of multiple (occa- 
sionally more than 4) bright, nearby companions. 

The final step in the determination of AGN+host mor- 
phologies involves comparing the Sersic-only and the 



Sersic+Point Source fits for each of the 57 sources in 
order to determine whether we significantly detect a cen- 
tral point source. Detailed simulations of the recovery of 
centr al point sources within ACS data (jSimmons fc Urrvl 
120081 ) suggest that failed detections of intrinsically- 
detectable point sources (i.e., false negatives) are rare, 
but spurious point-source detection (i.e., false positive) 
rates are strongly dependent on host morphology, rang- 
ing from ~ 1% in a sample of pure disk galaxies to ~ 25% 
for a sample of pure bulges. When GALFIT does not de- 
tect a central point source, therefore, this is most likely 
a robust indicator that any light from a central AGN is 
below the detection threshold of the image. Detection of 
a central point source in a Sersic+Point Source fit, on 
the other hand, is not a guarantee of its authenticity. 

We therefore assess the robustness of each point-source 
detection based on a number of criteria. First, if an F- 
test indicates that the addition of a point source sig- 
nificantly improves the xt goodness-of-fit parameter, 
we consider the point source genuine. Alternatively, 
we consider a point source genuine if the Sersic-only 
fit converges to an unrealistically high value for the 
Sersic index n (generally galaxies h ave n < 10; e.g., 
ICaon. Capaccioli. fe D'Onofrio|[T993l) . Lastly, we exam- 
ine the central portion of the fit residuals. We consider 
a point source genuine if the residual of the Sersic-only 
fit shows obvious signs of additional light from an un- 
resolved source (however, in practice at least one of the 
other two criteria also apply if this last criterion is true). 

Figure [1] shows a selection of fitted AGN+hosts; we 
present the full catalog of best-fit morphological pa- 
rameters for each source in Table [T] The uncertain- 
ties reported in Table [1] are determined by combining 
in quadrature the reported computational uncertainties 
from GALFIT and the uncertainties in morphological pa- 
rameter recovery determined f rom the simulation of ove r 
50,000 AGN host galaxies by iSimmons fc Urrvl (|200l) . 
The latter typically dominate the overall uncertainty in 
a morphological fit. 

3. RESULTS 

3.1. AGN Host Morphologies 

Just under half (25) of the host galaxies in this sam- 
ple are best described as single-component galaxies with 
smooth residuals after subtracting a single-component 
host fit. Fourteen have clumpy morphologies, with mul- 
tiple offset components either embedded in a larger com- 
ponent or very closely separated. Eighteen host galaxies 
(~ 32%) show signs of asymmetric features or tidal tails 
in the images and/or fit residuals; 8 of these appear to be 
involved in major mergers, consistent with an indepen- 
dent assessme nt of an overlapping sample using visual 



dent assessme nt ot an overlap ping s 
morphologies (jKocevski et al.ll2011[ ). 

Using the F-test to assess the significance of improve- 
ments to the fit with the addition of a central point 
source, we find that 90% of AGN+hosts have point 
source detections with at least 3er significance. The 
majority are at least 5a point-source detections. We 
quote point-source-subtracted host morphologies for the 
51 sources that have a point source at a minimum signif- 
icance level of 3er; for the remaining six sources we quote 
the morphology from the Sersic-only fit. As expected, 
for those sources where the addition of a point source 
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57 AGN Hosts: 
39 with 1.25 < 
28 with 1.76 < 



1.76 ■ 
2.67 




Fig. 2. — Distribution of fitted Sersic indices for the sample of 57 
AGN host galaxies at 1.25 < z < 2.67. The sample is characterized 
by disky morphologies (median n = 1.5). This is conservative: ex- 
tensiv e simulations of AGN host morphologies (Simmons & Urry 
120081) show that while a fitted n < 1.5 strictly confines the bulge 
contribution to less than 20%, AGN host galaxie s with fitted n = 4 
can still have significant intrinsic disks. We follow Simmons & Urry 
in separating host galaxies into three groups based on cuts at 
n = 1.5 and 3 (dotted lines). Sersic indices for sub-samples in 
two redshift bins spanning equal time intervals (At a 1.3 Gyr) are 
consistent with being drawn from the same parent sample. 

was not significant, the host morphology does not signif- 
icantly change between the Sersic only and the Sersic + 
point source fit . 

We show the fitted Sersic indices for the AGN host 
galaxies in Figure [2] Morphologies of the hosts, or (for 
the clumpy hosts) of the component within the host asso- 
ciated with the nuclear source, are typically disky. Just 
over half the sample has n < 1.5 (29 of 57 hosts), an un- 
ambiguously disk-dominated morphology. An additional 
15 have 1.5 < n < 3, which is consistent with interme- 
diate morphol ogies (bulge-to-total ratio s between 20 and 
65%). Indeed,! immons fc Urrvi ()2008f ) have shown that 
even galaxies fit with n = 4 can have a significant disk 
component (up to 45%). Five sources have 3 < n < 4, 
and 8 have n > 4. The overall sample can therefore be 
characterized as having a mix of morphological types, 
but disks are present in roughly three-quarters of the 
sample, and ~ 50% of hosts are consistent with having 
a bulge fraction of less than 20%. (In fact, 21 of the 
29 disk-dominated sources have n < 1, consistent with a 
bulge-to-total ratio of < 10%.) 

The redshift range over which the ^1251^ and F160W 
WFC3 filters sample the rest-frame B band spans ap- 
proximately 2.6 Gyr of cosmic time. Splitting the sample 
into two sub-samples of At « 1.3 Gyr each reveals no dif- 
ference in the distribution of morphologies between the 
two bins, according to a Kolmogorov-Smirnov (K-S) test. 

3.2. Host Galaxy Stellar Masses and Colors 

We determine stellar masses of the host galaxies using 
the multi-wavelength GOODS+CANDELS BVIzJHKs 
SEDs. We minimize the effect of the AGN on the host- 
galaxy SED fitting by first subtracting the point-source 
contribution to the flux in each HST band before fitting 
th e SEDs with the stellar populat ion synthesis tem plates 
of lMarastonl {2005) using FAST (IKriek et all 12009ft . 

Subtracting the point sources in the BVIzJH bands 
generally has only a small effect on the recovered galaxy 
masses and colors due to the fact that the host galax- 
ies dominate the observed SEDs: the median Lhost = 
10 x Lagn m JH . However, the effect is systematic: 
because using removing the AGN flux lowers the source 



flux used to calculate mass, the recovered stellar mass is 
lower than it would be were the presence of the AGN not 
accounted for. Additionally, for a sample such as this, 
where there is little bias against obscured AGN due to 
hard X-ray selection, but unobscured AGN are selected 
against by the requirement that the host galaxy mor- 
phology be recoverable, most nuclear point sources are 
reddened. Removing the AGN flux from the SED there- 
fore tends to recover a slightly bluer color than would be 
calculated if the effect of the AGN were neglected. This 
is likewise a small but systematic effect. 

Once the point-source SED is subtracted from the HST 
bands, extrapolating it to remove an estimated contribu- 
tion to the ifs-band flux makes no difference within the 
uncertainties. We therefore subtract the nuclear point- 
source contribution to the SED in the HST bands where 
we have enough angular resolution to confidently sepa- 
rate AGN from host galaxy, and refrain from unnecessary 
extrapolation of the SED to the ground-based Ks band. 

To test the dependence of recovered stellar masses on 
the choice of template, we c ompared stellar masses cal cu- 
lated using the templates of lBruzual fc Chariot ( 2003) to 
the masses calculated from iMarastonl (|2005| ) templates; 
we found that the recovered stellar mass is not strongly 
dependent on the choice of template set. The distribu- 
tion of derived masses is shown in Figure EH The min- 
imum and maximum masses of the sample are 2 x 10 8 
and 6 x 10 11 M Q , respectively; the average and median 
masses are 1.7 x 10 10 and 2.3 x 10 10 Mq, respectively. 
70% of the sample has M* > 10 10 M , which is ex- 
pected since AGN are t ypically found in massive hosts 
(|Cardamone et all 120101) : this probably reflects a selec- 
tion effect missing slow-growing black holes i n low-mass 
galax ies rather than an intrinsic distribution (lAird et al.l 
I2012T ). The average uncertainty in the mass is 0.8 dex. 

We visually inspected each of the best-fit templates 
from FAST and verified the goodness-of-fit. We use the 
best-fit templates to determine the rest -frame (U — V) 
colors of each host galaxy, follow i ng iBrammer et all 



((2009T) 



in assuming a 



Calzetti et all ([2000D dust law to 
correct for dust, such that (U — V) corr = (U — V) — 
AAy, where AAy = 0A7Ay- We show the color-mass 
diagrams for the AGN host sample in Figure 21 com- 
pared to inactive galaxies with the same redshifts (from 
the N EWFIRM medium-band survey; IWhitakc r et al.1 
I2011D . We separate the sample into two redshift bins 
divided at z = 1.76, so that each bin represents a time 
interval of At « 1.3 Gyr. 

A K-S test indicates no significant difference between 
the stellar mass distributions of the AGN hosts in the two 
redshift bins. The uncorrected U — V colors between the 
bins are inconsistent with being drawn from the same 
parent sample at the 97% level, but the significance is 
reduced to 75% (~ la) once the dust correction is ap- 
plied: correcting for dust shifts more galaxies into the 
blue cloud in the lower redshift bin than in the higher. 
A detailed comparison of the colors of morphological sub- 
samples in different redshift bins is highly uncertain due 
to small numbers in each sub-sample. However, we ob- 
serve intrinsically red disk-dominated hosts in each bin 
(at 85% confidence given the uncertainties), as well as 
blue strong-bulge hosts (99% confidence). 

Note that, although the stellar mass determination is 
unaffected (within the uncertainties) by the choice of 
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Fig. 3. — Histogram of stellar masses for AGN hosts with 
1.25 < z < 2.67. The full sample (solid black line) has a median 
value of 2 X 10 10 M (dotted line); we also show histograms for 
the unambiguously disk-dominated (n < 1.5; blue long-dash line), 
intermediate (1.5 < n < 3; green short-dash line), and strong- 
bulge (n > 3; red dotted line) sub-samples. Disk-dominated host 
galaxies tend to have lower masses than those with strong-bulge 
morphologies (K-S « 93%). 
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Fig. 4. — Rest-frame (U — V)ab color versus stellar mass for the 
host galaxies of AGN with 1.25 < z < 2.67, separated into two red- 
shift bins of approximately equal At 1.3 Gyr. Disk-dominated 
host galaxies (blue circles) populate both the "blue cloud" and 
"red sequence" portions of the diagram. Strong-bulge hosts (red 
squares) arc more likely to have redder colors, but not exclusively. 
The stellar mass distributions are statistically indistinguishable, 
and the dust-corrected color distributions of AGN host galaxies 
between the two redshift bins are only different at the 75% level, a 
marginal difference at best. Median error bars are shown in each 
panel. The size of each point representing a host galaxy is pro- 
portional to its AGN bolometric luminosity. Contours and gray 
point s show the positions o f inactive galaxies in the same redshift 
bins qWhitaker et alj|20Tl] ). 

model templates, the value of the dust correction Ay can 
vary significantly depending on the choice of template. 
We include differences in each s o urce's dust- correction 
betwe en the iBruzual fc Chariot] (|2003f ) and iMarastonl 
(|2005l) templates in the uncertainties shown in Figure 
21 noting that the choice of template does not change 
the overall result. 



3.3. Black Hole Masses and Accretion Rates 

Black hole mass, one of the fundamental properties of 
an AGN, can be measured directly only for very nearby 
sources. If a black hole is unobscured, measurements 
of broad spectral lines from the region near the black 
hole allow for calculations of central masses via the virial 
method. However, this is not generally possible for more 
obscured sources like those that dominate the present 
sample: only 1 AGN in this sample has a measurable 
broad line (see below), so this method is not feasible 
for the entire sample. Instead, we use well-characterized 
relations between host galaxy properties and black hole 
masses, which allow estimates of black hole masses out 
to high redshift as long as evolutionary corrections are 
considered. 

Locally, black hole ma ss correlates with th e ve- 
locity dispersi on (e.g., iFerrarese fc Merrittl 120001: 



Gebhardt etldl [2000). mass (e.g., iMagorrian et al.l 



19981: lHaring fc Rb3 12004ft and optical lumino s ity [e. H ., 
Marconi fc Huntl 120031 : IFerrarese fc Fordl 2005: Graham! 



a result of multiple effects (Treu et al. 2004: Borys et al. 


2005 


Woo et al.ll2005HPeng et al.ll2006t lAlexander et al. 


2008 


Woo et al.l 20081 Jahnke et al. 2009F Decarli et al. 


201C 


Merloni et all, 2010: Bluck et al.ll2011l). It arrears 



that bulge stellar populations are younger at earlier 
times, and it also appears that bulge growth may trail 
behind black hole growth at early times. However, the 
evoluti onary corrections are highly uncertain. R ecent 
work (jCisternas et al.l 120111: iMerloni et ail 120101 ) sug- 
gests that the relationship between the central black 
hole mass and the total stellar mass of the host galaxy 
does not evolve to z ~ 2. Other groups find that 
the relation evolves by a fact or of ~ 2 — 4 t o z ~ 2 
(jTrakhtenbrot fc Netzerl l2lHl iBennert et all 1201 lbj) . 
but the highest-redshift measurement to date of a 
single black hole with spatially-resolved spectroscopy 
(jlnskip et al.1 120 111 ) is consistent with the local black 
holc-galaxy relation. 

There are therefore many sources of uncertainty in the 
estimation of z ~ 2 black hole masses using galaxy-black 
hole correlations. In order to accurately characterize the 
uncertainties of our estimated black hole masses, we cal- 
culate black hole masses from host galaxy stellar masses 
via Monte Carlo simulations that account for the un- 
certainties in each quantity involved in the calculation 
— stellar mass, redshift and a possible evolutionary cor- 
rection of up to a factor of 3 — as well as the scatter 
in the black hole m ass-host stellar mass correlation of 
lHaring fc Rixl (|2004l ). We simulate 10 5 observations via 
Gaussian-random sampling within the uncertainties for 
each source and report the median black hole masses and 
la widths of the distributions of each source as our un- 
certainties. 

Table [1] gives the black hole mass estimates for each 
of our sources. We find that the black hole masses 
span a large range, from 2+\ 9 5 x 10 5 Af Q to 2^8 x 
10 9 M Q , with median and mean masses of 4 x 10 7 A/ Q 
and 3 x 10 7 M Q , respectively. Typical uncertainties are 
—0.7 and +1.0 dex, with asymmetries in the uncertain- 
ties primarily due to uncertainties in the host stellar 
mas ses, which are also calculated by Monte Carlo meth- 
ods (jKriek et al.ll2009T ). As a sanity check, we determined 
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57 AGN Hosts with 1.25 < z < 2.67 
• n < 1.5 
_ * 1.5 < n < 3 
■ n > 3 




10 4 



10 7 1C 
M BH [M ] 

Fig. 5. — Bolometric luminosity vs. black hole mass for AGN in 
disk- dominated (blue circles), intermediate (green triangles) and 
strong-bulge (red squares) host galaxies at 1.25 < z < 2.7 in this 
sample. Median error bars are shown at left beneath the legend. 
Lines shown indicate constant Eddington ratio: L = L Edd (solid 
line), L = OAL^dd (dashed line), and L = O.OlLgdd (dotted line). 
Growing black holes at these redshifts span three orders of magni- 
tude of growth rate, from Eddington-limited (and possibly super- 
Eddington) rates to highly sub-Eddington. One-third of the sample 
has L/L Ed d < 0.01; another ~ 30% have 0.01 < L/L EdA < 0.1. 
Bolometric luminosity is uncorrelated with host galaxy morphol- 
ogy- 

the black hole mass of X577 using the FWHM of the 
Mgll broad line and the c ontinuum luminosity at 3000 A 
(iMcLure fc Dunlodl200l. Th e latter is uncertain in the 
spectrum (|Szokolv et al.ll2004D . but the calculated mass 

of l.OtoJ x 1q8 M is fully consistent with that deter- 
mined using the Monte Carlo method described above. 

In addition to black hole masses, bolometric (total) lu- 
minosities are useful quantities because they allow for an 
estimate of the mass accretion ra te onto the black hole 
(modulo the radiative efficiency). iSimmons et al.l (|2011[ ) 
have shown that the model-dependent b olometric cor- 
rectio ns to the hard X-ray luminosity of iTreister et al.l 
are in very good agreement with bolometric lu- 
minosities recovered for obscured AGN by leveraging the 
observed point-source luminosities to extract the SEDs 
from the far-infrared to X-ray of each AGN. We there- 
fore use the absorption-corrected X-ray luminosity of 
each source, a long with the bolometric corrections of 
ITreister et al.1 (j2009f ). to estimate the bolometric lumi- 
nosities of each source. Given the X-ray luminosities in 
the sample, the bolometric correction ranges from 25 to 
76. 

With estimates of both bolometric luminosity and 
black hole mass, we explore the growth rates of the black 
holes in the sample in Figure [S] The growth rates of 
moderate-luminosity AGN at z ~ 2 span several or- 
ders of magnitude, from super-Eddington accretion to 
lQ~ 3 LEdd, with average and median L/LEdd values of 
0.06 and 0.05, respectively. One-quarter of the sam- 
ple (19 of 57 AGN) has L/L Edd < 0.01; 21 AGN have 



0.01 < L/LEdd < 0.1. The sample cannot therefore be 
characterized by a single growth ratc/Eddington ratio, 
but at least half of the sample is growing at less than 
10% of the Eddington (maximum) rate (96% confidence 
level) . 

4. DISCUSSION 

Because moderate-luminosity AGN far outnumber 
more luminous objects like quasars, they make up a 
large fraction of bla ck hole growth in the universe 
(jHasinger et al.l 120051 ) . They therefore represent the 
most common mode of black hole feeding that we can 
observe, and their feedback mechanism(s) represents the 
typical means of black hole-galaxy co-evolution. 

The host galaxies of the sample are generally charac- 
terized by strong disk contributions. Fully half of the 
sample of host galaxies is consistent with having less 
than 20% bulge contribution, and a further one-quarter 
of the AGN are hosted in galaxies with between 20 and 
65% of their light from a bulge. Very few host galax- 
ies in the sample are unamb iguously consisten t with 
the label "bulge-dominated": ISimmons fc Urrvl (|2008D 
showed that even a host galaxy with a fitted Sersic 
index of n > 4 (traditionally associated with a pure 
bulge or elliptical; there arc 8 in the sample) may have 
a strong disk contribution (up to 45% in some cases). 
Non-detection of a bulge, on the other hand, is a re- 
liable result according to multiple sets of host galaxy 
simulation s (Sanch ez et alJl2004t ISimmons fc UrrvH 2008: 
iGabor et aLll2009l : iPierce et al.ll2010aj) . 

That 75% of moderate-luminosity z ~ 2 AGN hosts in 
the CDF-S have strong disk c omponents agrees quanti- 



tatively with the pilot stud y of ISchawinski et al.l (|2011l) . 
Using visual classifications. IKocevski et al .1 (|2011l ) find a 
slightly lower fraction of disks (51%). However, visual 
classifications of host galaxies with even faint central nu- 
clear emission may be biased t oward recovering mor e 
bulge-like host morphologies (ISimmons et al.l I2012D . 
Given that 90% of moderate-luminosity AGN+hosts at 
z oj 2 have detected centra l point sources, the disk frac- 
tion of the IKocevski et al.l sample should be considered 
a lower limit. 

The high fraction of disky morphologies in moderate- 
luminosity AGN hosts at z ~ 2 (Figure [5]) means the 
majority of the host galaxies in the sample have assem- 
bled without building up prominent bulges. This im- 
plies that major mergers, which tend to destroy disks and 
lead to elliptical or bulge-dominated morphologies (e.g . , 
iToomrd [19771 : IWalker et al.1 [19961: IMartig et al.l 1201^ . 
play a small role in the growth of these galaxies and 
their central black holes. 

Note that, while a disk may be able to re-for m af- 
ter a gas-rich merger (e.g., iHopkins et all [2009), this 
merger process still strongly favors the formation of a 
bulge. Such a formation history is unlikely to be typical 
for the half of the sample with n < 1.5, and even less 
likely for the one-third of the sample with n < 1, which 
is consistent at these redshifts with a bulge fraction of 
less than 10% l|Simmons fc UrrvH2008l ) . Additionally, the 
majority of the sample does not show the signatures of a 
significant recent merger (the brighter features of such 
signatures could be d etected in images at this depth; 
ISchawinski et al.l I2012T ) . Thus gas-rich mergers that en- 
able disks to re-form yet fail to form significant bulges 
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either occurred at z ^ 2, so that by z ~ 2 they have 
relaxed to the point where they are not detected in deep 
rest-frame optical HST images, or they do not occur in 
moderate-luminosity AGN. 

Those galaxies that do appear to be highly disturbed 
also do not host the highest-luminosity AGN, as might 
be expected if the processes that disrupted the host 
also fed large amounts of material onto the black hole. 
Black hole-galaxy co-evolution scenarios in which ma- 
jor mergers trigger accretion onto the central super- 
massi ve black hole dSanders et al.lll988t iDi Matteo et al.l 
[200^lHoDk ins et al.ll2006al ). leading from an Eddington- 
limitcd quasar phase in a merger remnant host to a longer 
phase of slow black hole growth with a moderate- or low- 
luminosity AGN hosted by an elliptical galaxy, are in- 
consistent with the morphologies of moderate-luminosity 
AGN hosts at this early epoch. 

Disk-dominated hosts arc more likely to be observed as 
blue than bulges: although there are a substantial num- 
ber of intrinsically red disks in the sample, uncorrected 
colors of the sub-samples with n < 1.5 and n > 3 are 
different at the 99% level, according to a K-S test. How- 
ever, this difference disappears once the dust correction 
is applied, primarily because several of the strong-bulge 
hosts appear to be intrinsically blue, but dust-reddened. 
Blue early-type galaxies have also been seen for AGN 
hosts in other surveys at lower redshifts ( Sanchez" et al.l 
l2004HPierce et alJl2010btlCardamone et al.ll2010f ). 

However, the rest-frame U — V colors of z ~ 2 AGN 
hosts are not uniformly consistent with the "blue cloud" 
population at the same redshift, as expected from galaxy 
formation simulations where mergers tri gger black hole 
growth and a burst of star formation (Hopkin s et al.l 
120061* [Smnerville et al.l l2008f >. Even after a dust cor- 
rection is applied, ~ 50% of host galaxies are red, in 
both redshift bins. 

Merger scenarios tying the evolutionary stage of a post- 
merger to the decaying accretion rate of its central su- 
permassivc black hole also predict a correlation between 
AGN luminosity and host galaxy morphology. However, 
z ~ 2 AGN bolometric luminosity does not depend on 
host galaxy morphology: as seen in Figure black holes 
growing at or near the Eddington limit are found in disk- 
dominated, intermediate, and strong-bulge host galaxies, 
as are the slowest-growing black holes in the sample. 

Nor does the bolometric luminosity depend on redshift: 
black holes at all redshifts in the sample grow at all rates 
we are able to detect. In Figure |6] we show the bolometric 
luminosity versus redshift for the moderate-luminosity 
AGN in this sample. We also plot the flux limit of the 
sample, which defines a redshift-dependent lower limit to 
the bolometric luminosity. 

Figure [6] shows that the observed AGN populate ev- 
ery part of the parameter space that we do not select 
against. The range of AGN luminosities sampled implies 
a range of mass accretion rates spanning over two or- 
ders of magn itude. For an ass umed radiative efficiency 
of e = 0.15 (jElvis et al.l 120021 ). the black holes in the 
sample are growing at a rate of no more than ~ 2Mq 
per year for the highest-luminosity black hole. Lower 
accretion rates are more typical of the sample, with the 
average and median values being 0.03 and 0.02M Q yr -1 , 
respectively, and the lowest-luminosity AGN in the sam- 
ple accreting at 0.003M Q yr" 1 . 
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Fig. 6. — Bolometric luminosity vs. redshift for moderate- 
luminosity AGN in the GOODS-S/CA NDELS field. Gi ven the 
conversion of the X-ray flux limit of the IXue et all Q201H ) survey 
(dashed line) to a bolometric limit, our sources fill the observable 
area. The bolometric luminosity is converted into a mass accretion 
rate on t he right vertical axis, assuming a radiative efficiency of 
e = 0.15 JElvis et al.ll2003') . 

Such characteristically low mass accretion rates imply 
that slow- and moderate-growth periods of a typical su- 
permassive black hole do not require large, bulk flows of 
matter into the central region of the host galaxy, as one 
might expect from merger-triggered gas infall. Indeed, 
the median rate of 0.02M Q yr _1 is roughly consistent 
with expected stellar mass loss from a pa ssively evolv- 
ing stellar po pulation dCiotti et al.l 119911 : for a recent 
treatment, see ICiotti fc Ostrikerll2012D . ~ 

5. CONCLUSIONS 

We describe 57 AGN host galaxies from 1.25 < z < 
2.67 in the GO ODS-South field s elected from the 4 Ms X- 
ray catalogs of IXue et al.l (|20ll . Using HST WFC3/IR 
F125W and F160W images from CANDELS, we per- 
form parametric morphological fits in the rest-frame B 
band. We additionally calculate stellar masses, black 
hole masses, bolometric AGN luminosities, and rest- 
frame host galaxy (U — V) colors for each source. 

Our major results are as follows: 

• Host galaxy morphologies span a range from 
disk-dominated to bulge-dominated, but approxi- 
mately 50% of the sample is unambiguously disk- 
dominated, and another ~ 25% of the host sample 
have fitted Sersic indices consistent with strong (if 
not dominant) disk components. 

• Host morphologies are statistically indistinguish- 
able between two redshift bins of approximately 
equal time intervals (At m 1.3 Gyr; 1.25 < z < 
1.76 and 1.76 < z < 2.67). 

• Uncorrected host galaxy colors between the two 
redshift bins do appear different, with blue colors 
more likely at higher redshifts, but the significance 
of this difference is reduced (from rj 97% to w 74% 
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according to a K-S test) once a dust correction is 
applied. 

• We detect central point sources in the rest-frame B 
band to at least 3er in 90% of sources. The nuclei 
are typically faint compared to the host galaxy; the 
majority of detections are greater than 5a. 

• Central black hole masses, estimated from stellar 
masses, span a wide range (from 10 5 — 10 9 M Q ), 
with a median mass of 4 x 10 7 M Q . 

• AGN bolomctric luminosities range from 2 x 10 43 to 
2 x 10 46 erg s _1 ; this is limited only by the flux lim- 
its of the 4 Ms Chandra data and expected detec- 
tion numbers based on the AGN luminos ity func- 
tion at these redshifts (jCroom c t al. 2004]). Essen- 
tially, we find AGN accreting at all levels we are 
capable of detecting, with no redshift trend. 

• AGN bolometric luminosity is uncorrelated with 
host stellar mass, color or morphology, in contrast 
to the predictions of some models of merger-driven 
accretion. 

Moderate-luminosity AGN, which collectively com- 
prise a substantial fraction of total black hole growth 
at these redshifts, span several orders of magnitude in 
growth rate and are hosted by galaxies with a range of 
star formation rates, as indicated by a range of dust- 
corrected colors. While their host galaxies also span 
a range of morphologies, strong disks dominate. That 
fully half of growing black holes at these luminosities are 



hosted in galaxies with a (conservative) minimum of 80% 
of their optical light coming from an undisturbed disk 
likely indicates the importance of secular processes and 
minor mergers (that do not produce a significant bulge) 
in the co-evolution of black holes and galaxies. 
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a X-ray IDs from lXue et~aT1 poTlfl 

k X-ray luminosity from 0.5 — 8 keV 

c C: Clumpy; A: Asymmetric; M: Merger. 

i3~n 



For details, see Section 
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